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Abstract The hadronic equation of state for a neutron star is discussed with a 
particular emphasis on the symmetry energy. The resuhs of several mi- 
croscopic approaches are compared and also a new calculation in terms 
of the self-consistent Green function method is presented. In addition 
possible constraints on the symmetry energy coming from empirical in- 
formation on the neutron skin of finite nuclei are considered. 

1. Introduction 

In describing properties of neutron stars the equation of state (EoS) 
plays a crucial role as an input. During this workshop the effects of QCD 
degrees of freedom which become important at higher densities have been 
discussed in great detail. However, a quantitative understanding of the 
EoS also at lower densities is a prerequisite for the description of neutron 
stars properties. Conventional nuclear physics methods (i.e., in terms of 
hadronic degrees of freedom) are still being improved both conceptually 
as well as numerically, but all microscopic many-body approaches do 
require some sort of approximation scheme. 

The aim of this contribution is to review the present understanding of 
the hadronic EoS, and the nuclear symmetry energy (SE) in particular. 
In the past the latter quantity has been computed mostly in terms of the 
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Brueckner-Hartree-Fock (BHF) approach. In order to get an idea about 
the accuracy of the BHF result we present a recent calculation [1] in 
terms of the self-consistent Green function method, that can be consid- 
ered as a generalization of the BHF approach. The results are compared 
with those from various many-body approaches, such as variational and 
relativistic mean field approaches. In view of the large spread in the the- 
oretical predictions we also examine possible constraints on the nuclear 
SE that may be obtained from information from finite nuclei (such the 
neutron skin). 

2. Equation of state and symmetry energy 

The central quantity which determines neutron star properties is the 
EoS, at T = specified by the energy density e(p). For densities above, 
say, p = 0.1 fm~^ one assumes a charge neutral uniform matter consisting 
of protons, neutrons, electrons and muons; the conditions imposed are 
charge neutrality, Pp = Pe + Pfx-> and beta equilibrium, Pn = IJ'p + IJ-e with 
Me = /^M- 

Given the energy density e = eTv + + C/i the total pressure, P = 
Pn + Pe + P^, is obtained as P = P{e/(P) = J2 PifJ'i ~ with the 
chemical potentials given by pi = 

To illustrate the present status of the situation in Fig. 1 a compi- 
lation made by the Stony Brook group [2] of the relation between the 
pressure and density for a wide variety of EoSs (but still a tiny fraction 
of all available calculations) is shown. This figure shows that there is an 
appreciable spreading in the calculations of the pressure. Qualitatively 
one can distinguish three diff'erent classes. First the EoSs that corre- 
spond to self-bound systems (such as "strange quark matter") have the 
property that the pressure does not vanish at zero density. Then there 
are two globally parallel bands. It appears that most conventional non- 
relativistic approaches fall in the lower (softer) band and the covariant 
results in the upper (more rigid) band. It is worth noting that even at 
densities around and below saturation, po ~ 0.17 fm~^, the predictions 
vary appreciably (about a factor 4). At higher densities some models 
show a sudden change in the slope of the pressure. This behavior can 
be ascribed to the onset of new degrees of freedom, e.g., the appearance 
of hyperons, kaons and/or quarks which lead to a softening. 

It is well known that at lower densities the properties of the EoS are 
primarily determined by the SE [2]. The latter is defined in terms of a 
Taylor series expansion of the energy per particle for nuclear matter in 
terms of the asymmetry parameter a = {N — Z)/A (or equivalently the 
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Figure 1. Pressure versus density for a variety of EoS; taken from [2] where explar- 
nation of symbols is given. 

proton fraction x = Z/A), 

E{p, a) = E{p, 0) + S2{p)a^ + -S4(p)a^ + . . . , (1) 

where ^2 is the quadratic SE. It has been shown [3, 4] that the deviation 
from the parabohc law in Eq. (1), i.e., the term corresponding to S4, is 
quite small. 

Near the saturation density po the quadratic SE is expanded as 

^ , , ld'^E(p,a), Po , , AK . 

^^f") = 2 = «^ + ^(^ - po) + 3^(0 - po)' + ■ ■ ■ ■ (2) 

The quantity 04 corresponds to the SE at equilibrium density and the 
slope parameter pQ governs the density dependence. 
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As a result the pressure can be written as 
P{p,x) = p^^Ekr^ = p2[E'{p, 1/2) + S',ip){l - 2xf + ...]. (3) 

By using beta equilibrium in a neutron star, = — /Xp = — ^^^r^ ~ 
S'2(p)(l — 2x), and the result for the electron chemical potential, //g = 
one finds the proton fraction at saturation density, to 
be quite small, xq 0.04. Hence the pressure at saturation density can 
in good approximation be expressed in terms of (the density dependence 
of) the SE 

P{po) = po(l - 2xo){poS'2{po){l - 2xo) + S2ipo)xo) - ^0-52(^0). (4) 

3. How well do we know the symmetry energy? 

In the following the present status of calculations of the SE is reviewed; 
first the ones in which a microscopic nucleon-nucleon (NN) interaction 
is used, followed by the phenomenological mean field approaches. In 
section 4 we discuss possible constraints that can be obtained from em- 
pirical information. 

3.1 The SE in the BHF scheme 

In the Brueckner-Hartree-Fock (BHF) approximation, the Brueckner- 
Bethe-Goldstone (BBG) hole-line expansion is truncated at the two-hole- 
line level [5] . The short-range NN repulsion is treated by a rcsummation 
of the particle-particle ladder diagrams into an effective interaction or 
G-matrix. Self-consistency is required at the level of the BHF single- 
particle spectrum e^^^{k), 

e^HF^^^ = |1 + ^ Re < kk'\G{u; = 6^^^(fe) + e''''^ {k'))\kk' > . 

k'<kp 

(5) 

In the standard choice BHF the sclf-consistcncy requirement (5) is re- 
stricted to hole states {k < kp, the Fermi momentum) only, while the 
free spectrum is kept for particle states k > kp- The resulting gap in 
the s.p. spectrum at = fci? is avoided in the continuous-choice BHF 
(ccBHF), where Eq. (5) is used for both hole and particle states. The 
continuous choice for the s.p. spectrum is closer in spirit to many-body 
Green's function perturbation theory (see below). Moreover, recent re- 
sults indicate [6, 7] that the contribution of higher-order terms in the 
hole-line expansion is considerably smaller if the continuous choice is 
used. 
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Although the BHF approach has several shortcomings it provides a 
numerically simple and convenient scheme to provide insight in some 
aspects of the symmetry energy. 

Decomposition of SE. Some insight into the microscopic origin of 
the SE can be obtained by examining the separate contributions to the 
kinetic and potential energy [3]. 

The results of a ccBHF calculation [1] with the Reid93 interaction, 
including partial waves with J < 4 in the calculation of the G-matrix are 
presented in Fig. 2, where the SE is decomposed into various contribu- 
tions shown as a function of the density p. The kinetic energy contribu- 
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Figure 2 The SE S2 (full 
line) and the contributions 
to ^2 from the kinetic 
(dash-dotted line) and 
potential energy (dashed 
line), calculated within a 
ccBHF scheme and using 
the Rcid9.3 interaction. 
Also shown (dotted lines) 
are the T = and T = 1 
components of the poten- 
tial energy contribution. 



tion, iSkin, to the SE in BHF is given by the free Fermi-gas expression (it 
differs from the standard one, which is based upon the derivative rather 
than the finite difference) 

S'kin = ^kin,PNM " ^kin,SNM = Y^^^^^^^'^^ ~ ^~^) ' 
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and it determines to a large extent the density dependence of the SE. 
In Fig. 2 we also show the symmetry potential S'pot = 5'2 — Skin, which 
is much flatter, and the contributions to S'pot from both the isoscalar 
(T = 0) and isovector (T = 1) components of the interaction. Over 
the considered density range Spot is dominated by the positive T = 
part. The T = partial waves, containing the tensor force in the ^Si-^Di 
channel which gives a major contribution to the potential energy in SNM, 
do not contribute to the PNM energy. The T = contribution peaks at 
p !=a 0.3 fm~^. The decrease of this contribution at higher densities is 
compensated by the increase of the T = 1 potential energy, with as a net 
result a much weaker density-dependence of the total potential energy. 

Dependence of the SE on NN interaction. Engvik et al. [8] 
have performed lowest-order BHF calculations in SNM and PNM for all 
"modern" potentials (CD-Bonn, Argonne vl8, Reid93, Nijmegen I and 
II), which fit the Nijmegen NN scattering database with high accuracy. 

They concluded that for small and normal densities the SE is largely 
independent of the interaction used, e.g. at po the value of 04 varies 
around an average value of 04=29.8 MeV by about 1 MeV. At larger 
densities the spread becomes larger; however, the SE keeps increasing 
with density, in contrast to some of the older potentials like Argonne 
vl4 and the original Reid interaction (Reid68) for which S{p) tends to 
saturate at densities larger than p = 0.4 fm~^. 

3.2 Variational approach 

Detailed studies for SNM and PNM using variational chain summa- 
tion (VCS) techniques were performed by Wiringa et al. [9] for the Ar- 
gonne Avl4 NN interaction in combination with the Urbana UVIII three- 
nucleon interaction (3NI), and by Akmal et al. [10] for the modern Avl8 
NN potential in combination with the UIX-3NI. 

It should be noted that the results of VCS and BHF calculations 
using the same NN interaction disagree in several aspects. For instance 
for SNM the VCS and BHF calculations saturate at different values of 
density [5]. As for the SE using only two-body interactions the VCS 
approach yields a smaller value for 04 than BHF (see Table 1), and as 
a function of density in VCS the SE levels off at p ~ 0.6 fm~^, whereas 
the BHF result continues to increase. Therefore it seems natural to ask 
whether the inclusion of more correlations by extending the BHF method 
(which is basically a mean field approximation) will lead to results closer 
to those of VCS. 
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3.3 Self-consistent Green function method 

As noted above BHF has several deficiencies: it does not saturate at 
the empirical density, and it violates the Hugenholz-van Hove theorem. 

In recent years several groups have considered the replacement of the 
BBG hole-line expansion with self-consistent Green's function (SCGF) 
theory [11-13]. In ref.[13, 14] the binding energy for SNM was calcu- 
lated within the SCGF framework using the Reid93 potential. In ref.[l] 
we have extended these calculations to PNM and considered the corre- 
sponding SE. 

The SCGF approach differs in two important ways from the BHF 
scheme. Firstly, within SCGF particles and holes are treated on an equal 
footing, whereas in BHF only intermediate particle (k > kp) states are 
included in the ladder diagrams. This aspect ensures thermodynamic 
consistency, e.g. the Fermi energy or chemical potential of the nucleons 
equals the binding energy at saturation (i.e. it fulfills the Hugenholz- 
van Hove theorem). In the low-density limit BHF and SCGF coincide. 
As the density increases the phase space for hole-hole propagation is no 
longer negligible, and this leads to an important repulsive effect on the 
total energy. Second, the SCGF generates realistic spectral functions, 
which are used to evaluate the effective interaction and the correspond- 
ing nucleon self-energy. The spectral functions exhibits a depletion of 
the quasi-particle peak and the appearance of single-particle strength at 
large values of energy and momentum, in agreement with experimental 
information from (e, e'p) reactions. This is in contrast with the BHF 
approach where all s.p. strength is concentrated at the BHF-energy as 
determined from Eq. (5). 

In the SCGF approach the particle states (k > kp), which are absent 
in the BHF energy sum rule do contribute according to the energy sum 
rule {d = 1(2) for PNM(SNM)) 

expressed in terms of the nucleon spectral function Sh{k,ui). 

To illustrate the difference between the ccBHF and SCGF approaches 
the results for both SNM and PNM are compared in the left and central 
panels of Fig. 3 for the Reid93 interaction. One sees that the inclusion 
of high-momentum nucleons leads roughly to a doubling of the kinetic 
and potential energy in SNM, as compared to BHF. The net effect for 
the total energy of SNM is a repulsion, increasing with density [13]. This 
leads to a stiffer equation of state, and a shift of the SNM saturation 
density towards lower densities. The above effects which are dominated 
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by the tensor force (the isoscalar ^Si-^Di partial wave) are much smaher 
in PNM. 

The corresponding SE, shown in the right panel of Fig. 3, is dominated 
by the shift in the total energy for SNM, and lies below the ccBHF 
symmetry energy in the entire density-range. At po = 0-16 fm~^ the 
parameter 04 is reduced from 28.9 McV to 24.9 MeV, while the slope pQ 
remains almost the same (2.0 MeVfm"^ compared to 1.9 MeVfm"^ in 
BHF). 
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Figure 3. The energy per particle for syrnrnctric nuclear matter (left panel) and 
pure neutron matter (central panel) for the Rcid93 interaction. The dashed line refers 
to a ccBHF calculation, the full line to a SCGF calculation. The right panel displays 
the SE in these two approaches. 



3.4 Three-body force in microscopic approaches 

As noted above at higher densities the EoS is sensitive to 3NF contri- 
butions. Whereas the 3NF for low densities seems now well understood 
its contribution to nuclear matter densities remains unsettled. In prac- 
tice in calculations of the symmetry energy in the BHF approach two 
types of 3NF have been used in calculations; in rcf.[4] the microscopic 
3NF based upon meson exchange by Grange et al. was used, and in ref. 
[15] as well in most VCS calculations the Urbana interaction. The latter 
has in addition to an attractive microscopic two-pion exchange part a 
repulsive phenomenological part constructed in such a way that the em- 
pirical saturation point for SNM is reproduced. Also in practice in the 
BHF approach to simplify the computational efforts the 3NF is reduced 
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to a density dependent two-body force by averaging over the position of 
the third particle. 

In general inclusion of the 3NF stiffens the SE for p > pQ. For example, 
with the AvlS+UIX interaction 04 increases from 27 to 30 MeV and po 
from 2 to 3 MeVfm-^ (see Table 1). 

In view of the important contribution from the 3NF to the SE at 
higher densities there seems a clear need to get a more quantitative 
understanding of the repulsive part of the 3NF in PNM. 

3.5 Dirac-Brueckner-Hartree-Fock 

The SE has also been computed in the Dirac-Brueckner-Hartree-Fock 
(DBHF) approach [16, 17]. Although some of the shortcomings of BHF 
(mean field treatment of the Pauli operator and violation of Hugenholz- 
van Hove theorem) persist in the DBHF method it appears that in the 
latter approach one is able to reproduce the saturation properties of SNM 
without the inclusion of a three-body interaction. A general feature of 
relativistic methods is that the SE increases almost linearly with density, 
and more rapidly than in the non-relativistic case. This difference can 
be attributed to two effects. First the covariant kinetic energy which is 

inversely proportional to \Jk'p + m*^ is larger because of the decreasing 
Dirac mass, m*, with increasing density. Secondly the contribution from 
rho-exchange appears to be larger than in the non-relativistic case [16]. 

3.6 Relativistic mean-field approach 

Relativistic mean field (RMF) models have been applied succesfully 
to describe properties of finite nuclei. In general ground state energies, 
spin-orbit splittings, etc. can be described well in terms of a few pa- 
rameters ref. [18]. Recently it has lead to the suggestion that the bulk 
SE is strongly correlated with the neutron skin [19, 20] (see below). In 
essence the method is based upon the use of energy-density functional 
(EDF) theory. 

In practice covariant approaches are formulated either in terms of 
a covariant lagrangian with a, co and p exchange (and possibly other 
mesons) [18, 21], or in terms of contact interactions [20], solved as an 
EDF in the Hartree-Fock approximation. Sets of model parameters are 
determined by fitting bound state properties of nuclei. Specifically the 
isovector degree of freedom is governed by the exchange of isovector 
mesons; in case of p-meson exchange the (positive definite) contribution 
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to the SE is given by 



04 = 



(8) 



and its potential energy contribution to which scales with that for 04, 

2 

is [20]. Typical values obtained for po are around 4-6 MeV fm~^, 
and 04 ~ 30-36 MeV, i.e., considerably larger than in non-relativistic 
approaches (a large part of the enhancement can be ascribed to the fact 
that the kinetic contribution is larger, because m* < ra). 

Recently this approach was extended by inclusion of the isovector- 
scalar partner, the 5- meson, of the isoscalar scalar cr- meson [22]. Unfor- 
tunately the value of the coupling for the (5-meson cannot be determined 
well by fitting properties of stable nuclei. Also in its simplest, density in- 
dependent form, the inclusion of the (5-meson leads to an even larger net 
value for po- This happens because of the presence of the Lorentz factor 
m* /E in the scalar potential contribution, ~ which decreases 

with increasing density. 



Recently the density dependence of the symmetry energy has been 
computed in chiral perturbation effective field theory, described by pions 
plus one cutoff parameter, A, to simulate the short distance behavior 
[23] . The nuclear matter calculations have been performed up to three- 
loop order; the density dependence comes from the replacement of the 
free nucleon propagator by the in-medium one, specified by the Fermi 
momentum kp 



The resulting EoS is expressed as an expansion in powers oi kp and 
the value of A 0.65 GeV is adjusted to the empirical binding energy 
per nucleon. In its present form the validity of this approach is clearly 
confined to relatively small values of the Fermi momentum, i.e. rather 
low densities. Remarkably for SNM the calculation appears to be able 
to reproduce the microscopic EoS up to p ~ 0.5 fm~^. As for the SE 
the value obtained in this approach for 04 = 33 MeV is in reasonable 
agreement with the empirical one; however, at higher densities [p > 0.2 
fm^^) a downward bending is predicted (see Fig. 4) which is not present 
in other approaches. 
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Effective field theory 
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3.8 Comparison of results 

To summarize the present situation in Fig. 4 the resulting density 
dependence of the SE for the approaches discussed above are compared 
(excluding the 3NF contribution). One sees that the covariant mod- 
els predict a much larger increase of the SE with the density than the 
non-relativistic approaches. The lowest-order BHF method predicts a 
somewhat higher value for 04 than both the VCS and SCGF methods, 
which lead to very similar results; whether that can be ascribed to a 
consistent treatment of correlations in these methods, or is fortuitous, is 
not clear. 
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Figure 4- Overview of several theoretical predictions for the SE: Brueckner-Hartree- 
Fock (continuous choice) with Reid93 potential (circles), self-consistent Green func- 
tion theory with Reid93 potential (full line), variational calculation from [9] with 
Argonne Avl4 potential (dashed line), DBHF calculation from [16] (triangles), rel- 
ativistic mean-field model from [22] (squares), effective field theory from [23] (dash- 
dotted line). 

Clearly at least part of the differences should be attributed to the 
different selection of the mesonic degrees of freedom in the various mod- 
els. In the microscopic approaches the tensor force mediated by vr and 
p exchange seems to play the dominant role. On the other hand in 
the mean field approach explicit pion exchange is usually not included 
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Table 1. Results for the symmetry energy parameters 04 (in MeV) and po (in 

MeVfm-3) 





BHF 


BHF+3NF[5] 


VCS[10] 


VCS+3NF[10] 


SCGF[1] 


DBHF[16] 


RMF[20] 




28.9 


30 


27 


30 


24.9 


31 


30-36 


Po 


1.9 


2.9 


2.0 


3.1 


2.0 


3.0 
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and hence there the isovector effect solely comes from the shorter range 

p-cxchange. In fact it has been argued that in contrast to isoscalar 
properties the long-range pion exchange should play an essential role in 
determining the isovector properties [20]. 

4. Empirical information on the SE 

In view of the existing uncertainties in the calculation of the SE one 
may ask whether from finite nuclei one can obtain experimental con- 
straints on the symmetry energy as a function of density. In this section 
some recent activities pertaining to this issue are reviewed. 

4.1 Relation between SE and neutron skin 

Recently in applying the non-relativistic Skyrme Hartree-Fock (SHF) 
model Brown [19] noted that certain combinations of parameters in the 
SHF are not well determined by a fit to ground state binding energies 
alone; as a result a wide range of predictions for the EoS for PNM can be 
obtained. At the same time he found a correlation between the derivative 
of the neutron star EoS (i.e., basically the symmetry pressure po) and 
the neutron skin in ^*^^Pb. 

Subsequently Furnstahl [20] in a more extensive study pointed out 
that within the framework of mean field models (both non-relativistic 
Skyrme as well as relativistic models) there exists an almost linear em- 
pirical correlation between theoretical predictions for both 04 and its 
density dependence, poi and the neutron skin, /S.R = Rn — Rp, in heavy 
nuclei. This is illustrated for ^°^Pb in Fig. 5 (from ref.[20]; a similar cor- 
relation is found between AR and po). Note that whereas the Skyrme 
results cover a wide range of AR values the RMF predictions in general 
lead to AR > 0.20 fm. 

The interpretation of the nuclear matter results in part depends on 
the question whether there is a surface contribution to the SE in finite 
nuclei. In ref.[24] it was found that for heavy nuclei the latter is of minor 
importance, which has also been confirmed in ref.[20]. 
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Figure 5. Neutron skin thickness versus 04 for Pb for a variety of mean field 
models (from [20]). The circles correspond to results for the Skyrme force, the squares 

to RMF models with mesons, and the triangles to RMF with point couplings; the 
shaded area indicates the range of A7i values consistent with the present empirical 
information for ^"^Pb. 



4.2 Insight in the correlation of AiZ and 04 

The above observation suggests an intriguing relationship between a 
bulk property of infinite nuclear matter and a surface property of finite 
systems. 

Here we want to point out that this correlation can be understood natu- 
rally in terms of the Landau- Migdal approach. To this end we consider a 

simple mean- field model (see, e.g., rcf. [25]) with the Hamiltonian consist- 
ing of the single-particle mean field part Hq and the residual particle-hole 
interaction Hph'- 

H = Y,{Ta + Uo{Xa) + Ui{Xa) + Uc{Xa)) + Hph, (9) 
a 

where (x = (r, a, r)) 

Uo{x) = Uo{r)+Uso{xy, Ui{x) = ^-Spot^r^^); Uc{x) = ic/c(r)(l-r(-^)). 

(10) 

Here, the mean field potential includes the phenomenological isoscalar 
part Uo{x) along with the isovector Ui{x) and the Coulomb Uc{x) 
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parts calculated consistently in the Hartree approximation; Uo{r) and 
Uso{x) = Uso{r)a ■ I are the central and spin-orbit parts of the isoscalar 
mean field, respectively, and S'pot('') is the potential part of the symme- 
try energy. 

In the Landau-Migdal approach the effective isovector particle-hole in- 
teraction Hpfi is given by 

Hph = + G'aaffb)ramra - n), (H) 

a>b 

where F' and G' are the phenomenological Landau-Migdal parameters. 

The model Hamiltonian H in Eq. (9) preserves isospin symmetry if the 
condition 

[H,f(-^] = U^j\ (12) 
is fulfilled, where f(-) = Eri~\ ^ = E^c(ra)ri~\ With the use 

a a 

of Eqs. (9), (11) the condition eq.(12) in the random phase approxima- 
tion (RPA) formalism leads to a self-consistency relation between the 
symmetry potential and the Landau parameter F' [26]: 

5pot(r) = 2F'n(-)(r), (13) 

where n^^\r) = n^{r) — nP{r) is the neutron excess density. Thus, in this 
model the depth of the symmetry potential is controlled by the Landau- 
Migdal parameter F' (analogously to the role played by the parameter 
gj^ in rclativistic mean field models). 

Spoti^) is obtained from Eq. (13) by an iterative procedure; the re- 
sulting dependence of Ai? on the dimensionless parameter /' = F' / (300 
MeVfm^) shown in fig. 6 indeed illustrates that Ai? depends almost lin- 
early on /'. Then with the use of the Migdal relation [27] which relates 
the SE and /', 

a4 = ^(1 + 2/0, (14) 

an almost linear correlation between the symmetry energy, 04, and the 
neutron skin is found. 

To get more insight in the role of /' we consider small variations 
5F' . Neglecting the variation of n'^^)(r) with respect to 5F' one has 
a linear variation of the symmetry potential: 5Spot{r) = 25F'n^~\r). 
Then in first order perturbation theory, such a variation of Spot causes 
the following variation of the ground-state wave function 
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Figure 6. Neutron skin in ^"^Pb vs. the Landau-Migdal parameter /'. 



with s labeling the eigenstates of the nuclear Hamiltonian and a single- 
particle operator iV^") = X] n*^~)(ra)Ti^'*. Consequently the variation of 

a 

the expectation value of the single-particle operator V^~^ = Y^r 

a 

with (0|F(-)|0) = NRl - ZR^ can be written as 

Re(0|iV(-)|s)(s|y(-)|0) 



R,6iAR)=6F'^Yl 



En — Es 



(16) 



In practice the sum in Eq. (16) is exhausted mainly by the isovec- 
tor monopole resonance of which the high excitation energy (about 24 
MeV in ^''^Pb) justifies the perturbative consideration. We checked that 
Eq. (16) is able to reproduce directly calculated d{AR) shown in fig. 6 
with the accuracy of about 10%. 

As for the correlation between AR and po one would need information 
on the density dependence of F'. Sofar as we know it has not been ex- 
tracted from data on stable nuclei. In the approximation that F' is 
density independent one naturally finds the po is proportional to 04. 

4.3 Experimental information on AR for ^°*Pb 

What are the experimental constraints on the neutron skin? A variety 
of experimental approaches have been explored in the past to obtain 
information on AR. To a certain extent all analysis contain a certain 
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model dependence, which is difficult to estimate quantitatively. It is 
not our intention to present a full overview of existing methods for the 
special case of ^'^^Pb. In particular the results obtained in the past from 
the analysis of elastic scattering of protons and neutrons have varied 
depending upon specifics of the analysis employed. At present the most 
accurate value for AR comes from a recent detailed analysis of the elastic 
proton scattering reaction at E = 0.5 — 1 GeV [28], and of neutron and 
proton scattering at £^ = 40 — 100 MeV [29]. For details we refer to 
these papers. Here we restrict ourselves to a discussion of the some less 
well known methods that have the potential to provide more accurate 
information on the neutron skin in the future. 

Anti-protonic atoms. Recently neutron density distributions in a 
series of nuclei were deduced from anti-protonic atoms [30]. The ba- 
sic method determines the ratio of neutron and proton distributions at 
large differences by means of a measurement of the annihilation prod- 
ucts which indicates whether the antiproton was captured on a neutron 
or a proton. In the analysis two assumptions are made. First a best fit 
value for the ratio Rj of the imaginary parts of the free space pp and 
pn scattering lengths equal to unity is adopted. Secondly in order to 
reduce the density ratio at the annihilation side to a a ratio of rms radii 
a two-parameter Fermi distribution is assumed. The model dependence 
introduced by these assumptions is difficult to judge. Since a large num- 
ber of nuclei have been measured one may argue that the value of Ri is 
fixed empirically. 

Parity violating electron scattering. Recently it has been pro- 
posed to use the (parity violating) weak interaction to probe the neutron 
distribution. This is probably the least model dependent approach [31]. 
The weak potential between electron and a nucleus 

y(r) = F(r)+75^(r), (17) 

where the axial potential A(r) = ■^f^pwi'i^)- The weak charge is mainly 
determined by neutrons 

pwir) = (1 -Asm'^9w)pp{r) - pn{r), (18) 

with sin^ 9w ~ 0.23. In a scattering experiment using polarised electrons 
one can determine the cross section asymmetry [31] which comes from 
the interference between the A and V contributions. Using the measured 
neutron form factor at small finite value of and the existing infor- 
mation on the charge distribution one can uniquely extract the neutron 
skin. Some slight model dependence comes from the need to assume a 
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certain radial dependence for the neutron density, to extract i?„ from a 
finite form factor. 

Giant dipole resonance. Isovector giant resonances contain in- 
formation about the SE through the restoring force. In particular the 
excitation of the isovector giant dipole resonance (GDR) with isoscalar 
probes has been used to extract /S.R/R [32]. In the distorted wave Born 
approximation optical model analysis of the cross section the neutron 
and proton transition densities are needed as an input. For example, in 
the Goldhaber- Teller picture these are expressed as 

with K the oscillation amplitude and {i = p,n); one assumes ground 
state neutron and proton distributions of the form (a; = (A'' — Z)/A) 

Pi{r) = 1(1 ± X T 7x)pir - c(l ± 7x/3)), (20) 

where 7 is related to the neutron skin AR, 7 = 2(iv-Z) ^-^/-^0" Then 
ioT N > Z the isovector transition density takes on the form 

A /X N — Z ,dp c d?p. 

In practice one studies the excitation of the GDR by alpha parti- 
cle scattering (isoscalar probe). By comparing the experimental cross 
section with the theoretical one (calculated as a function of the ratio 
AR/R) the value of AR can be deduced [32]. 

It is difficult to make a quantitative estimate of the uncertainty in 
the result coming from the model dependence of the approach. In the 
analysis several assumptions must be made, such as the radial shape 
of the density oscillations and the actual values of the optical model 
parameters. 

We note that also other types of isovector giant resonances have been 
suggested as a source of information on the neutron skin, such as the 
spin-dipole giant resonance [33] and the isobaric analogue state [34] . At 
present studies of these reactions have not led to quantitative constraints 
for the neutron skin of ^'^^Pb. 



Results for A.R. In table 2 we present a summary of some recent 
results on AR in ^''^Pb. One sees that the recent results are consistent 
with AR values in the range 0.10-0.16 fm. It appears from Fig. 5 that 
this range is consistent with the conventional Skyrme model approach 
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but tends to disagree with the results of the RMF models considered in 
[20]. Also from the correlation plot between AR and po shown in Fig. 11 
in ref.[20] one may conclude that a small value for po ~ 2.0 MeVfm"^ is 
favoured over the larger values from covariant models. 



Table 2. Summary of recent results for AR in ^°*Pb 



method 


AT? [fm] 


error [fm] 


ref 


{p,p') at 0.5-1.04 GeV 


0.097 


0.014 


[28] 


nucleon scattering (40-200 MeV) 


0.17 




[29] 


anti-protonic atoms 


0.15 


0.02 


[30] 


giant dipole resonance excitation 


0.19 


0.09 


[32] 


parity violating electron scattering 


planned 


1% 


[31] 



4.4 Information on the SE from heavy-ion 
reactions 

In principle the density dependence of the SE at higher densities (and 
further away from N = Z) can be probed by means of heavy-ion reacties 
using neutron rich radioactive beams. In ref. [35] possible observable 
effects from the isovector field are considered in terms of the RMF model. 
Of particular interest seems the study of the ratio of tt~^ /tt~ cross sections 
and its energy dependence which appears sensitive to details of the SE 
[36]. 

5. Constraints on EoS from neutron stars 

Given the EoS as a function of density the mass versus radius relation 
of a neutron star can be obtained in the standard way with the use 
of the Tolman-Oppenheimer-Volkov equation. Do we need the EoS up 
to all densities or can we already draw some conclusions from the lower 
densities? In [2] it was argued that there exists a quantitative correlation 
between the neutron star radius and the pressure which does not depend 
strongly on the EoS at the highest densities. Generally speaking a stiff 
(soft) EoS implies a large (small) radius. 

Clearly a simultaneous measurement of the radius and the mass of 
the same star would discriminate between various EoS. The presently 
available observations do not yet offer strong constraints. For example, 
most observed star masses fall in the range of 1.3-1.5 Mq. These values 
seem to rule out the very soft EoSs predicted with hyperons present, 
which yield M^ax ~ 1.3M0 [5]. This may suggest that either the NY 
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and YY interactions used as an input must be reexamined, or that the 
use of the BHF approach at high densities is not reliable. 

Other constraints come from recent observations from X-ray satellites. 
Most robust seem the data from the low mass X-ray binary EXO 0478- 
676 obtained by Cottam et al. [37] . Prom the redshifted absorption lines 
from ionized Fe and O a gravitational redshift z = 0.23 was deduced; 
this gives rise to a mass-to-radius relation 

M/Mo = (1 - (21) 

with RgQ = 2.95 km. As discussed in detail in various contributions 
in this workshop at present this constraint appears fully consistent not 
only with conventional EoSs, but also with most more exotic ones. 
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